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Abstract 

We propose and study a local map capable of describing the full variety of dynamical states, 
ranging from regular to chaotic, obtained when a nematic liquid crystal is subjected to a steady 
shear flow. The map is formulated in terms of a quaternion parametrization of rotations of the local 
frame described by the axes of the nematic director, subdirector and the joint normal to these, with 
two additional scalars describing the strength of ordering. Our model yields kayaking, wagging, 
tumbling, aligned and coexistence states, in agreement with previous formulations based on coupled 
ordinary differential equations. Such a map can serve as a building block for the construction of 
lattice models of the complex spatio-temporal states predicted for sheared nematics. 
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Driven complex fluids exhibit an unusual variety of dynamical statesQ, [2], [3, Q, [|], [f3, Q] . 
When such fluids are sheared uniformly, the stress response is regular at very small shear 
rates. However, at larger shear rates the response is often intrinsically unsteady, exhibiting 
oscillations in space and time as a prelude to intermittency and chaos^, Q, Such chaos 
associated with rheological response or "rheochaos" , occurs in regimes where the Reynolds 
number is very small. It must thus be a consequence of constitutive and not convective 
non-linearities, originating in the coupling of the flow to structural or orientational variables 
describing the local state of the fluid |9j, [lOj. The diverse possibilities for internal degrees of 
freedom in complex fluids, such as the orientation of nematogenic molecules, layer stacking 
in lamellar and onion phases and heterogeneities arising from local jamming in colloidal 
suspensions, implies that the study of the rheology of complex fluids should illuminate a 
variety of non-trivial steady states in driven soft matter. 

Recent rheological studies of "living polymers" , solutions of worm-like micelles in which 



the energies for scission and recombination are thermally accessible, obtain an oscr 
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latory 

Q. it 



response to steady shear at low shear rates which turns chaotic at larger shear rates 
has been argued that a hydrodynamic description of this behaviour requires a field describing 
the local orientation of the polymer, motivating a treatment of the problem of an orientable 
fluid, such nematic, in a uniform shear flowl 
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, Q, [l^ . Nonlinear relaxation equations 



for the symmetric, traceless second rank tensor Q characterizing local order in a sheared 

15, Q, 0, Q. Assuming spatial uniformity, a 



nematic have been derived 
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system of 5 coupled ordinary differential equations (ODEs) for the 5 independent components 
of Q in a suitable tensor basis is obtained. Solving this system of equations yields a complex 
phase diagram admitting many states - aligned, tumbling, wagging, kayak-wagging, kayak- 
tumbling and chaotic - as functions of the shear rate 7 and aphenomenological relaxation 



time which is a parameter in the equations of motion 
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211 ] . Recent work adds spatial 



variations: numerical studies of the partial differential equations thus obtained vie 



diagram containing spatio-temporally regular, intermittent and chaotic states 22 



d a phase 



23]. 



The degrees of freedom which enter a coarse-grained description of an orientable fluid are 
mesoscopic. Spatio-temporal structure arises from the coupling of locally ordered regions, 
through processes such as molecular diffusion, flow-induced dissipation and advection. A 
powerful approach to understanding complex spatio-temporal dynamics is based on the 
study of coupled map lattices, a numerical scheme in which maps placed on the sites of a 
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lattice evolve both via local dynamics as well as through couplings to neighbouring sites [241]. 
However, the utility of this methodology in a specific context is often severely limited by the 
availability of local maps able to describe the spatially uniform case. This paper addresses 
this requirement in the context of a model for rheochaos, proposing the first local map 
description of the regular and chaotic states obtained in sheared nematics. 

There is, in general, no systematic procedure for the construction of such maps. However, 
it is reasonable to require that any such map should accurately reproduce the full variety of 
states obtained through the study of the corresponding ODEs. It should also enable useful 
physical insights through a sensible choice of physical variables. One obvious possibility is 
simply the discretization of the governing ODEs. Such a choice of variables, however, is not 
particularly illuminating as these equations are formulated in terms of the components of Q 
in a specific space-fixed tensor basis, rather than in terms of variables more natural to the 
problem. 

We have thus explored an alternative formulation of this problem, constructing a local 
map in terms of quaternion variables. These variables encode the dynamics of the orthog- 
onal set of axes associated with the eigenvectors of Q, i.e. the director, sub-director and 
the joint normal to these. Our approach incorporates biaxiality, is formulated in terms of 
physically accessible variables and is computationally straightforward to implement. Our 
results, summarized in the phase diagram of Fig. [I], are in good agreement with previous 



work based on ODEs [2l|], but provide an efficient alternative to such methods 25] 
Defining b := |(b + b T ) — ^(trh)5 to be the symmetric-traceless 



part o: 



tensor b, the equation of motion for Q in a passive velocity field is [U 



2l|: 



the second-rank 



^ - 2fTQ - 2af^Q + t ^ = -V2^T (1) 

at v r a 

where the tensor O = |((Vv) T — Vv), T = |((Vv) T + Vv) and Vv is the velocity gradient 
tensor, with v = 7?/e x , where e x is a unit vector in the x— direction. The velocity is along the 
x direction, the velocity gradient is along the y direction, while z is the vorticity direction. 
The quantities r a > and r ap are phenomenological quantities related to relaxation times, 
a describes the change of alignment caused by T and 3> = <90/<9Q, with 0(Q) = ^AQ : Q — 
|\/6-B(Q ■ Q) : Q+|C(Q : Q) 2 . The notation Q : Q represents QijQji, with repeated indices 
summed over. Here A = A Q (1 — T*/T), and B and C are constrained by the conditions 
A > 0, B > 0, C > and B 2 > \A C. Scaling t = t*r a /A k , v = v*A k /r a and a = a*a k , 
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FIG. 1: [Color online] The phase diagram of steady states in our model, illustrating regimes in 
which the following steady states are obtained for a generic initial condition: an aligned state 
denoted tumbling state labelled as 'T', a wagging state 'W, a kayak-tumbling state 'KT', 

a kayak-wagging state denoted by 'KW and a complex state denoted as 'C. This phase diagram 
closely resembles phase diagrams plotted in Refs. [211 ] . 



Eqn. (P) can be written in dimensionless form, ^-2n^*-2af\Q* + (eQ*-3VQQ*^Q* + 
2(Q* : Q*)Q*) = yj\ \ k T* where A k = A Q (1 - T*/T k ) = 2B 2 /9C,a k = a eq {T k ) = 2B/3C is 
the (nonzero) equilibrium value of the scalar order parameter a at the transition temperature 
T k , \ k = — fv^T^: and 9 = (1 — ^r)/ (1 — jr) is the reduced temperature. 

The Q tensor admits the following parametrization: Qij = ^ \niUj — + 
Y (jriimj — lilj), where s% and S2 represent the magnitude of the ordering along n (the 
director) and m (the subdirector), with n and m unit vectors and 1 = n x m. The dynam- 
ics of Q thus involves both the dynamics of the frame defined by n, m and 1 as well as the 
dynamics of s\ and S2- The frame dynamics can be represented in many equivalent ways, 
such as through coordinate matrices, axis-angle or Euler angle representations. However, 
the coordinate matrix representation requires a large number of parameters, the axis-angle 
representation suffers from redundancy and the use of the Euler-angle representation is 
marred by the "gimbal-lock" problem 26j . Our parametrization of the frame dynamics uses 
quaternion variables, providing an elegant, compact and numerically stable alternative to 
these representations. 

Equations for n, m and 1 as well as for the order parameter amplitudes s\ and s'2 can be 
derived by considering a reference frame in which the director and subdirector are stationary 
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FIG. 2: [Color Online] The sequence of three main panels shows the power spectrum associated with 
states in the regimes labelled (a) T and (b) W in the phase diagram of Fig. [TJ The topmost panel 
(c) shows a mixed state (M) (not shown separately in Fig. [I]), associated with the boundary between 
W and T The inset labelled (i) in all these panels shows typical plots of the time-dependence of 
the z-component of the director n z and the angle <p made by the projection of the director on the 
x — y plane with the x— axis. The insets labelled (ii) in all these panels show the trajectory in the 
s\ — S2 plane. 



(body frame). In the body frame, denoted by primed vectors, the director can be chosen to 
be n' = (1, 0, 0), the subdirector to be m' = (0, 1, 0), with 1' = (0, 0, 1). The transformation 
matrix A which maps vectors from the lab frame to the body frame, can be denned in 
terms of quaternion parameters (eo, . . . , e^j constrained by e\ + e\ + e\ + e| = 1[27|. The 
quantities n = (n x , n y , n z ), m = (m x ,m y ,m z ) and 1 = (l x ,l y ,l z ) are easily obtained using 
this mapping, yielding ODE's for the parameters Si, s 2 , e , ei, e 2 , e 3 , e 4 . These are converted 
into a map using a first-order Euler scheme. After each discrete time step, we renormalise 
the quaternion variable. Choosing a and 9 equal to zero for all the results reported here in 
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common with earlier work, our map is then defined through 
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We choose A = 0.01 for all our calculations. (The phase boundaries shown in Fig. [1] exhibit 
a weak dependence on At. However, provided At is chosen small enough, this dependence 
may be neglected.) The superscript 't' indicates that the values of the variables are taken 
at the t'th discrete time step. The control parameters are the dimensionless shear rate 7 
and Afc. In place of the 5 coupled ODE's used in the conventional parametrization of the 
dynamics of Q, we have 6 equations constrained by the normalization requirement, thereby 
conserving the number of degrees of freedom. 

In our numerical analysis of the map, we start typically from random initial conditions, 
omitting sufficient transients (~ 10 5 time steps) to ensure that the asymptotic attractor of 
the dynamics is reached. Our analysis includes inspection of the (i) power spectrum, (ii) 
phase portraits, (iii) bifurcation diagrams and (iv) time series of the different relevant vari- 
ables. Figs. [2] and [3] show the variety of states obtained in our numerical calculations. Each 
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FIG. 3: [Color online] The sequence of three main panels shows the power spectrum associated 
with states in the regimes labelled (a) KT (kayak-tumbling), (b) KW (kayak- wagging) and (c) C 
(complex or chaotic) in the phase diagram of Fig. [TJ The inset labelled (i) in all these panels shows 
typical plots of the time-dependence of the z-component of the director n z and the angle 4> made 
by the projection of the director on the x — y plane with the x— axis. The insets labelled (ii) in 
all these panels show the trajectory in the s\ — S2 plane. 

sub- figure, labelled as Figs. [2] (a) - (c) and Figs. [3] (a)-(c), has the following structure: The 
first inset, labelled (i) for all figures, describes the time dependence of n z , the z-component 
of the director, and the angle <fi made by the projection of the director on the x—y plane with 
the x— axis. The second inset, labelled (ii) for all figures, plots the quantities measuring 
the amount of ordering along director and sub-director against each other, providing the 
attractor of the system in the Si — S2 plane for a generic initial condition. The main panel 
in each of the sub-figures shows the power spectrum of s%, ln(|v4(/)| 2 ) against frequency / 
on a semi-log plot. 

The following states are easily identified: (I) An Aligned state denoted as 'A' in the 
phase-diagram of Fig. [TJ but omitted, for brevity, from the states shown in Fig. [2] and Fig. [3j 
In the aligned state, neither the frame orientation, nor s± and S2, vary in time. The director 
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is aligned with the flow at a fixed angle; (II) A Tumbling state, in which the director lies 
in the shear plane (the xy plane) and rotates about the vorticity direction (the z axis). 
Fig. [2(a) (i) indicates that this state is a stable in-plane state, since the z-component of the 
director is zero. Also, the angle made by the projection of the director on the x-y plane 
varies smoothly between 7r/2 and -n/2. Fig. 12(a) (ii) shows the periodic character of this 
state. This state is labelled as 'T' in the phase-diagram of Fig{TJ (III) A Wagging state, in 
which the director lies in the shear plane, but oscillates between two values. Note that Fig. 
[2j(b)(i) indicates that this state is a stable in-plane state. Also, the director oscillates back 
and forth in-plane as indicated in Fig. [2j(b)(ii). Fig. [2j(b) shows that this state is a periodic 
state with sharp delta-function peaks in the power spectrum. These states are denoted as 
'W in the phase-diagram in Fig{TJ In addition, we obtain (IV) A Kayak- Tumbling state, 
equivalent to the tumbling state, but in which the director is out of the shear plane. Thus, as 
shown in Fig. [3(a) n z ^ and the projection of the director on the xy plane rotates through 
a full cycle. Such states are temporally periodic, as shown in Fig. [3(a); the regular cycles 
evident in the map of s± vs. S2 (Fig. [3(a) (ii)) is a further indication of periodic behaviour. 
These states are noted as 'KT' in the phase-diagram of Fig. [TJ (V) A Kayak- Wagging 
state where, as in KT, the director is out of plane, but the projection of the director on the 
shear plane oscillates between two values. The properties of such states are illustrated in 
Fig. [3(b). Such states are again temporally periodic. The cyclic trajectory of the system in 
the Si — S2 plane (Fig. [3(b) (ii)) further confirms such periodic behaviour. These states are 
denoted by 'KW' in the phase-diagram of Fig. CH (VI) A Mixed state, typically found close 
to the boundaries between wagging and tumbling states, whose properties are illustrated 
in Fig. [2(c) . In such states, the director exhibits both oscillation and complete rotations. 
Power spectra obtained at the boundaries of this phase, for example near Xk = 0.99 and 
7 = 4.0, have a broad range of frequencies, and, (VII) A Complex state, in which the 
director lies out of the shear plane but both oscillates and rotates. The complex phase 
exhibits chaotic behaviour, as can be seen in Fig. [3(c). Note that the delta function peaks 
in the power spectrum exhibited by the periodic states discussed earlier have broadened into 
a continuum of frequencies. The plot of Si vs. S2 displays no regular structure. These state 
are noted as 'C in the phase-diagram in Fig{TJ In addition to these states, we also obtain a 
log rolling state in which the director is perpendicular to the shear plane (not shown). 
The range of dynamical states manifest in this problem is also evident from bifurcation 
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diagrams obtained at fixed 7 = 4.0, varying Such a cut intersects KT, T, W, KW, C and 
A states in the phase diagram. The quantities n z and the Poincare section of s\ show that 
n z = for the T, W and A states, while the KT, KW and C states are out-of-plane states 
with n z 7^ 0. Further, the S\ section, displays a fixed point for the aligned state, regular 
cycles for the KT, T W and KW states and irregular (chaotic) behavior for the C state. 

Aradian and Cates have recently studied a minimal model for rheochaos in shear- 
thickening fluids, using equations which describe a shear-banding system coupled to a 
retarded stress response 29j. These authors connect their model system to a modified 
Fitzhugh-Nagumo model, a dynamical system with a variety of interesting and complex 
phases. Fielding and Olmsted study instabilities in shear-thinning fluids, where the instabil- 
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ity originates in the multi-branched character of the constitutive relation 
al. report a study of the PDE's describing the dynamics of Q, characterizing spatio-temporal 
routes to chaotic behaviour in sheared nematics 22j. These studies allow for spatial variation 



- although restricted so far to the one- dimensional case - whereas our local map describes the 
spatially uniform situation. However, the dynamical system we study is obtained directly 



3,Q. 



Whether cou- 



from the underlying dynamics, in contrast to the approaches of Refs. 
pling maps of the sort we construct permits a complete description of the spatio-temporal 
structure obtained in Ref. 22| remains to be seen. 



In conclusion, we have proposed a local map describing the variety of dynamical states 
obtained in a model for sheared nematics. Our phase diagram, Fig. [H contains all non-trivial 
dynamical states obtained in previous work. It also closely resembles, even quantitatively, 
phase diagrams obtained in previous work which used ordinary differential equations formu- 
lated in continuous time. Our work thus supplies a crucial ingredient in the construction of 
coupled map lattice approaches to the spatio-temporal aspects of rheological chaos, a prob- 
lem currently at the boundaries of our understanding of the dynamics of complex fluids. 
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